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For a long time it was known that deuteron, as a weakly coupled nucleon pair, has no excited states. 
However, A.M. Baldin et al, commenting results of the first physical experiment with accelerated 
nuclei at JINR synchrophasotron, assumed as far back as in 1979 that one of peaks in a differential 
cross-section may arise due to an ’’excited state of deuterium”. We have established that one of 
the peaks in the cross-section may be explained indeed in this way and corresponds to the dibaryon 
reported by WASA-at-COSY Collaboration. Another peak in the same region is interpreted most 
likely by interference of several ^‘-resonances, and this possibility was also mentioned in the paper 
by A.M. Baldin et al. Further experimental studies based on modern experimental facilities and 
more abundant statistics are necessary to verify these observations. 

PACS numbers: 25.45.De, 25.10.+S, 27.10.+h 


I. INTRODUCTION 

Recently a proposal of QCD investigation at high den¬ 
sity and low temperature, complementary to the heavy 
nuclear collisions, was suggestedi*^. The proposal is 
based on the fact that a large number of nucleons in the 
interaction region is not necessary for the phase transi¬ 
tion to occur, and only a change of the vacuum state 
should be initiated by some experimental environment. 
In particular, observation of multi-baryons (MB) may be 
a direct evidence of phase transitions in small nucleon 
systems. Separation of a MB mass from the secondary 
particle background is feasible if the MB decay width is 
narrow enough. That requires the excitation energy of 
MB produced should be low. For this purpose, it is rea¬ 
sonable to select only those experimental events in which 
the MB creation is accompanied with a high momentum 
particle, taking away an essential part of the energy from 
the interaction region - method of cumulative particle^. 
In this paper, we focus on a verification of this concept 
by the use of older experimental data taken at JINR syn¬ 
chrophasotron. 

An experiment^ was designed for measurement of 
cross-sections of pp-, ND-, and DD-interactions at 8.9 
GeV momentum of primary protons and deuterons. A 
model of the detector operation was briefly described in^. 
Its parameters were established by means of measuring 
differential elastic cross-sections for proton-proton scat¬ 
tering in a known kinematic region. Three peaks were 
observed in the spectrum of the missing masses of the 
reaction D-fD—J-Mj^-I-D at t = —0.495 GeV^. Here we 
shall concentrate only on one of them, called the third 
peak in the original paper. In regard to the third peak, 
M.A. Baldin et al suggested that it might occur due to: 
a contribution of an excited state of deuteron; scatter¬ 
ing of a constituent quark (entering into the composi¬ 
tion of the incident deuteron) by target deuteron; and, 


in addition, N*-baryon production. Experimental find¬ 
ings occurred after the paper— was written give a cause 
for re-examination of the suggestions mentioned above. 

The present paper might be considered as a particular 
proposal for experimental search of phase transitions in 
small nucleon systems. 

II. CONSTITUENT QUARK SCATTERING 

Elastic scattering of a constituent quark by the target 
deuteron may be considered in the framework of a model 
in which values of momentum and mass of the projectile 
quark are considered in the form 

Pq = xPi, Mq = xMD, 

where x is determined from kinematics of the reaction. 
A necessary relation between quark mass and known ex¬ 
perimental parameters can be found as follows. Let us 
denote by 1-1-2 3-1-4 a reaction at issue, where the pro¬ 

jectile, target and registered particles are designated by 
1, 2 and 4, correspondingly, and 3 denotes an object X 
which mass should be determined. Two different expres¬ 
sions for the Lorentz invariant Mandelstam variable u, 
u = {pi- piY and u = {p 2 - Ps)^, where pi = 
allow to connect Mx and cos 9 which describs the escape 
direction of the particle 4 in the laboratory system. A 
value energy of particle 4 as function of M 2 , M 4 and t 
may be found by making use of a relation t = (p 2 — paY- 
In addition, = Ei+ E 2 — E 4 . Proceeding on this way, 
one obtains 


Eit-\- |Pi I t(— 4M|, -I- t) cos 9 ’ 

and Mq = 0.351 GeV for cos0 = 0.396. This number 
contradicts manifestly to estimations of modern quark 
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TABLE I: Spin, parity and width of N* included in our PWA. 
The data are given by Particle Data Group^. 


N* 

Sn* 

Pn* 

Tn*, MeV 

N(1440) 

1/2 

1 

300 

N(1520) 

3/2 

-1 

115 

N(1535) 

1/2 

-1 

150 


models: see, e.g.;^ where Mq = 0.318 GeV. On the other 
hand, we shall see below that a peak at cos0 = 0.396 
corresponds remarkably to the dibaryon found by WASA- 
at-COSY Collaboration^. 


III. PARTIAL-WAVE ANALYSIS (PWA) AND 
SU{6) (g) 0(3) QUARK SPECTROSCOPY 

Now let us turn to study of a possible contribution of 
N + D ^ N* + D reactions to the experimental cross- 
section. Isotopic spin conservation constrains isospin of 
N* to be equal to 1/2. Therefore, A—baryon excitations 
of nucleon may be ignored here, and among N* excita¬ 
tions only N(1440), N(1520) and N(1535) are important 
in the kinematic region under consideration. Main char¬ 
acteristics of the baryon resonances taken into account 
are shown in Table HI 

Besides the spatial parity, P, conservation, one should 
respect the angular momentum, J = L-f S, preservation. 
In sp-approximation, appropriate to the hard collisions, 
only L = 0 and L = 1 eigenvalues of the orbital momen¬ 
tum can be considered. In these terms, parities of initial 
and final states may be expressed as follows: 

P^ = PnPd{-1)^^ = (- 1 )^^ =Pf = ( 1 ) 

Further PWA may be simplified essentially via appli¬ 
cation of the 51/(6) 0 0(3) description of baryon exci¬ 
tations, suggested by R.H. Dalitz and co-authorsi. Ac¬ 
cording to it, spin Sjv* of a nucleon resonance N* may 
be represented as follows: 

S^. = + 1, (2) 

where Sjv is spin of the unexcited nucleon N and 1 is 
orbital momentum of quarks inside of the excited nu¬ 
cleon N*. Using o, it is readily seen that for each 
partial wave, which is characterized by fixed values of 
J and P, a value of parity Pn* of nucleon resonance N* 
determines totally a possible behavior of I and L val¬ 
ues. For N(1440), one has Pyr* = Pn = 1 which implies 
I = 0, and, subject to ([T|), also L/ = Li. For N(1520) 
and N(1535), Pyr* = — 1; therefore 1 = 1. According to 
(ED and Table m we can interpret spins of N(1520) and 
N(1535) as two different manners of summation, using 
Clebsch-Gordan coefficients, of quark orbital momentum 
I = 1 and initial spin Sn=1/‘^ of unexcited nucleon. The 


parity conservation leads to simultaneous change of L 
and I values in two possible ways: 

Lj = 1 —^ Lf =0, /i = 0 —>■ 1/ = 1, (3) 

and 

Pi = 0 —>■ Lf = 1, li = 0 ^ If = 1. (4) 

In the frame of SU{ 6 ) 0 0(3) spectroscopy, these cases 
correspond to conservation of eigenvalues of operator 

= (L -I- 1)^, which are equal to 2 and 0, accordingly. 
Operator of the total orbital momentum M commutes 
with M^, and we can develop a more detail picture in¬ 
cluding account of a direction of M. Below we shall 
consider centrally symmetric interaction conserving the 
direction of M. In this case, conservation of the total 
angular and orbital momenta implies preservation of the 
total spin, S = J — M, of the system and our description 
admits a further development. 

A general expression of the N + D ^ N* + D ampli¬ 
tude linear relative to ^at, Sd and invariant under time 
reversal and space rotation or reflection is as follows^ 

T(SN,SD) = Ci+C 2 (SN + SD)-j + C3(SN-SD)-i (5) 

+C4(SN-j)(^D-i) + C5(SN-'k)(SD-k) + Ce(SN-i){SD-i), 

where 

._pxp' , p-p' . p + p' 

•’“IpxpT “Ip-pT '~|p + pr 

p and p' are momenta of the ingoing nucleon and outgo¬ 
ing N*. Here Ci are scalar functions which may depend 
only on a scalar (p-p 0 /|p||p 1 which is in one-to-one cor¬ 
respondence with cos 9. In fact, we should claim C 3 = 0, 
for Sn — Su does not commute with (Sjv + S^)^ and 
the corresponding term breaks conservation of an abso¬ 
lute value of the total spin. Similarly, it is possible to 
show that C 4 = C 5 = Cq = Because of the total 
spin conservation, a term proportional to (Sat -I- S^)^ is 
not included in as far as it is proportional to unit 
operator for any state with a total spin S fixed (where, 
as usual, = S{S + 1)). Efficiently, it is included in Ci. 

Thus, we have seen that the parity conservation ad¬ 
mits concordant alteration of L and I according to © 
and ©. From the physical point of view (|3]) corresponds 
to swapping external orbital momentum of N-l-D system 
into nucleon, and corresponds to excitation of the 
both external, L, and intranucleonic, 1, momenta. These 
processes may be described by a nonlocal operator (R-r) 
included in the interaction amplitude. Here R is a polar 
vector given in the laboratory system, which is directed 
at center of inertia of N-l-D system, and r is a polar vec¬ 
tor pointed at center of mass of the nucleon colliding with 
deuteron. Without loss of generality, we may also sug¬ 
gest (R-R) = (r-r) = 1. Then T-matrix describing 
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production of baryon from Table |T] may be written in the 
form 

T{N + D^N* + D)=A + B{SN + SD)-i 

+ {-R-r)[C + D{SN + SD)-j]. (6) 

Here A and B describe spin independent and spin de¬ 
pendent parts of interaction corresponding to N(1440) 
production. Similarly, C and D describe interaction cor¬ 
responding to N(1520) and N(1535). Using an identity 

(R • r) = i {R+r- + R-r+) + R^r^, 

and well-known formulae for R and r operators^ 


^ f{S, 3/2) (C + DjSr, + S^) ■ j) r(jV(1520) ^ N + n) 

-^.^(1520) “ ~ *-^Af(1520)rAr(i520) 

( 7 ) 

^ f{S, 1/2) {C + D{Sn + Sz,) ■ j) r(fV(1535) ^ N + ir) 

^N{1535) ~ ~ *-^Af(1535)rAf(1535) 

Analogous expressions take place for N-l-D—^N-|-7r-|-7r-|-D 
and N-l-D—^-N-I- t^-I-D transitions. In ([7]), scalar functions 
A, B, C, D are the same as in and factors f{S, N*) 
may be found on basis of Clebsch-Gordan coefficients, 
as it was mentioned in previous section. Following this 
prescription, one can find 

ns,s^,)= E E E 

CTi—±1/2 o'2—0,±1 m—0,±1 


(L = 1, M = 0| |L = 0, M = 0) = -i/\/3, 
{I = l,m = 0|rz \l = 0,m = 0) = —ifVS, 


1 


-tTllCT2 


S,ai + a2) { -cTilm 


SN*,cri + m 


and 


{L=1,M = \L = 0,M = 0) = 



2 + V2, f 
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(Z = 1, m = -|-1| r+ |Z = 0, m = 0) = +i\/2j3, 

it is possible to find that amplitudes of the processes dH) 
and O are equal to 1 and 1/3, accordingly. 


IV. OBSERVABLE PARTICLES, 
CROSS-SECTION 

In fact, baryon resonances N(1440), N(1520) and 
N(1535) were not observed directly. They were present 
in an intermediate state and may be identified only via 
their decay products. Therefore interference terms cor¬ 
responding simultaneous propagation of matter through 
several quantum states with different spins and pari¬ 
ties should be taken into account. We take for granted 
that possible final states tolerating macroscopic recog¬ 
nition may contain Ntt, Ntttt and Nr;, of course, besides 
deuteron. For N(1440) and N(1520), corresponding decay 
probabilities can be estimated as wit^ ~ 0.65, W 2 tv ~ 0.35, 
Wff ~ 0; and witt ~ 0.5, W 2 tt ~ 0.1, Wrj ~ 0.4 for N(1535), 
see^. 

Baryon resonances leave imprint of their existence only 
as propagators in total amplitude. For example, a tran¬ 
sition N-l-D—^N-I- tt-I-D is described by the following T- 
matrix: 


T{N + D ^ N + TT + D) = 




4(72 + ^3 + 76), 


where we adopted notations of the Clebsch-Gordan coef¬ 
ficients from^. 

Here we should re-arrange a usual formula for cross- 
section^, 

(Pa IT 1 

dt dM]^ Ai/2(s, vPj) {2 Sn + 1)(2S'd + 1) 


X E] / dLips(M|-, decay product s)\TMiMi\^ , (8) 

Mi.Mj 

where and M f are spin projections of particles in ini¬ 
tial and final states, into terms of our model of the orbital 
nucleon excitations. To this end, we replace averaging 
over Mi and summation over Mj by corresponding op¬ 
eration over Si and E/, which are total spin projections 
of quarks in initial and final states. For nonpolarized ini¬ 
tial states, probabilities of occurrence of ^ = 1/2 and 
5 = 3/2 are equal to 1/3 and 2/3, accordingly. Taking 
into account that the contribution of orbital excitations 
is already included by means of /(5, Sn*), we may write: 

(25^ + l)(25z, + l) '' 


(A + H(Sjv + Sd) ■ j)T(V(1440) ^ V + tt) 
-^.^(1440) ~ ^X — tM7v(l440)r7V(1440) 


= 3 E |7SiE/P + 

S/=±l/2 


3 E ITe.E/P: 

i:/=±l/2,±3/2 
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and then transform 

EPW = E = E 

= (^E,E. = ^^Tr(Trt). 

Now it is easy to prove a relation 

(25^ + 1)(25b + 1) '' 


= ^ E Tr(r(5)T(5)t), 

C_ 1 3 
2 ’ 2 

which means that values of total spin 5" = 1/2 and 
3/2, as well as all its projections E = ±1/2 and E = 
±1/2, ±3/2, correspondingly, give equal contribution to 
the final result. It should be stressed that the sign 1 of 
Hermitian conjugation refers to T as to spin operator, 
and it does not mean transposition of other variables^!. 

Calculations of Tr(TTl) may be completed with mak¬ 
ing use of relations: Tr(S • j) = 0, 


Tr(l) 


2, 5=1/2, 
4, 5 = 3/2, 


Tr(S.j)2 


1/2, 5=1/2, 
5, 5 = 3/2. 


Absolute values of the decay amplitudes are fixed in 
terms of decay widths^, 

= 2ii^ / ^Lips(Mi, /) ^ \T{N* ^ f)\\ 

'' Mf 


where subscript N* denotes a particular baryon reso¬ 
nance, / is its decay products. We confine our estima¬ 
tions of interference between different baryon resonances 
to operations with phase space averaged values. For this 
purpose, we define^^ 

2M^t ^ /)PLips(Mi, /), 

and substituted^ 

( 25 ^. ± i)(|T(iv/^/)p) (|r(iv; ^ /)p) 


xe*(“’““^)Lips(M|:,/) = (9) 


for 


/ dLips(Mi, /) ^ T{N: ^ f)T*iN; ^ /) 

if Mx is greater than N* decay threshold and zero oth¬ 
erwise. Here baryon resonances are different, i ^ and 


decay particles are the same for the both multipliers un¬ 
der integral sign. 

Strictly speaking, separate control of spin projections 
of N* is not kept in mind in our description, but only 
projection of total spin of quarks in the final state of re¬ 
action N±D—>'N*±D. Therefore we should take into ac¬ 
count availability of deuteron too and replace Mj with 
E/ and Sx* with 5 in the previous formulae. Such a 
treatment may be understood as summation over quark 
spin projections inside N* and spectator deuteron. Con¬ 
tribution of orbital excitations into spin projection of N* 
is already included explicitly by means of /(5, 5Ar*), as 
it was mentioned above. This new interpretation of spin 
summation rule is an inevitable corollary of considera¬ 
tion of baryon as a compound system with its own inner 
structure. 

In the accepted approximation, only phases of the de¬ 
cay amplitudes cXi may be used as adjustable parameters 
for experimental data matching. In addition, eight real 
numbers corresponding complex parameters A, B, C, D 
in 0 are brought into play for this purpose. Interference 
terms corresponding decays of N* via r] are absent since 
cross-sections of this channel are negligible quantities but 
for one of the resonances under consideration (see values 
Wri in beginning of this section). The final formula de¬ 
scribing the experimental data may be written in the 
following form: 


(Pa TT 

dt dM^ 


/ dLips(Mi,/)Tr(T(5,/)Tt(5,/))±A, 

f,s 

where / = Ntt, Nr], Ninr, 5 = |, | and additional ad¬ 
justable parameter E describes a contribution of direct 
pion production near = 1.5 2 GeV^. 


V. SOME DETAILS OF NUMERICAL 
CALCULATIONS 

To reach an optimum in describing the experimental 
data we minimized total deviation square for 22 exper¬ 
imental points from the theoretical curve. Ten central 
experimental points, as the most important, were taken 
with unit weights and six ones on their left and six ones 
on their right were scaled with 0.5 significance. MAPLE 
procedure NLPSolve for the local minimum search was 
used for optimal selection of theoretical parameters. Sev¬ 
eral series, each containing 20 000 different sets of random 
initial values of parameters, were generated and only 30 
percent of them were finished without interruption be¬ 
cause of very big number of steps towards a local mini¬ 
mum. Points of the interruptions were considered as local 
minima too, because they usually correspond to wander¬ 
ings along valleys. Then the best local minimum was 
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taken for each of the series, and values of objective func¬ 
tion corresponding to them were compared. They turned 
out to be equal within accuracy of 11 decimal digits. All 
the best local optima have demonstrated that experimen¬ 
tal data demand unambiguously: 

|A|=0, |C|=0. (10) 

This means that phases (j)A and (f>c of complex numbers 
A and C have no impact upon objective function. For 
removal of degeneration, we have fixed (pA = (l^c = 0 and 
introduced condition (nni) explicitly into minimizing func¬ 
tional. Now the normal mode of NLPSolve performance 
increased up to 55 percent signalling, nevertheless, that a 
large degeneration still persisted. Three series of numeri¬ 
cal experiments, containing 100, 1000 and 20 000 events, 
with random selections of initial values of the remaining 
parameters were fulfilled. They showed that parameters 
|i?|, In I and E are identical in all the cases and are deter¬ 
mined with accuracy of 4 and 6 decimal digits already in 
the series with 100 and 1000 events. However, all phases 
underwent rather strong changes with growth of statis¬ 
tics, signalling that minimizing functional remains degen¬ 
erate with respect to them. Thus, the optimization prob¬ 
lem does not allow us to determine phases of parameters 
A, B, C, D, TiN* N+tt) and T(7V; ^ N+tt+tt), be¬ 
cause many of their sets describe equally well the exper¬ 
imental data. A grade of fidelity of reproduction of the 
experimental data by this model may be seen in Fig. 1. 

We have also fulfilled evaluation of the model param¬ 
eters using only 10 experimental points taken straight 
from the fine structure location, trying to enhance an 
impact of the most important region. It was technically 
fully regular procedure, as far as we had only 8 indepen¬ 
dent parameters at that stage. However, an agreement 
between theory and experiment has not been improved 
even in this case. 


VI. CONCLUSIONS AND DISCUSSION 

Numerical analysis fulfilled within the bounds of our 
model has revealed two nonobvious properties of hard 
N-D and D-D scattering. First of all, it was estab¬ 
lished that experimental datai^ show strong spin depen¬ 
dence of N+D—>'N*-|-D transition amplitude, see ([6]) with 
A = C = 0. Secondly, comparison of the experimental 
data and theory shown in Fig. 1 makes an explicit hint 
of dibaryon production in this kinematic region. Indeed, 
on the one hand, consideration only usual nucleon ex¬ 
citations cannot explain the fine structure shown in the 
figure. On the other hand, assumption about presence 
of a dibaryon at M 2 B ~ 2.38 GeV, F 2 B ~ 70 MeV, 
seen by WASA-at-COSY Collaboration^ allows one to 
explain it very naturally. Isospin conservation predicts 
certainly that reaction D-l-D —dibaryon -|- D should 
yield dibaryon with isospin / = 0, which also corresponds 
to the WASA-at-COSY result^. Thus, our consideration 
of the data on the hard deuteron-deuteron scattering^ 



FIG. 1: The experimental data (bars) in the range of the 
third peak and their explanation by the sum of contribu¬ 
tions of N-l-D—>-N*-|-D reactions (dashed line). The lower scale 
corresponds to the kinematics of reaction N-l-D—>-N*-|-D, the 
top one describes reactions D-l-D—>-X-|-D, which implies the 
dibaryon production. A possible contribution of a dibaryon 
at 2.37 GeV, T rs 70 MeV reported by WASA-at-GOSY 
Gollaboratioii^ into the cross-section is shown with the over¬ 
turned solid line. 

meets the expectation to observe the transition of nucleon 
matter into other states using the method of cumulative 
particle which allows to recognize quasi-resonance peaks 
in the reaction cross-section. 

To check our conclusions, it would be enough to mea¬ 
sure with a good precision production cross-sections of 
N(1440), N(1520) and N(1535) from N-hD ^ N*-hD reac¬ 
tions in appropriate kinematic region, and direct produc¬ 
tion of pions therein. This allows one to take into account 
the background. In addition, repeating experiment^ with 
higher accuracy is necessary too for unambiguous recog¬ 
nition of dibaryon by its mass and width. Theoretical 
and experimental study of the phases entering into ex¬ 
pression for production amplitude is ineffectual in this 
respect, so long as resultant cross-section is weakly de¬ 
pendent on them (see previous section). Investigation 
of decay products of dibaryon will make it possible to 
identify its spin and parity and compare with = d"*" 
observed in^. 

It is interesting to review ability of the lattice QCD 
to say something definite about existence of dibaryons. 
All lattice QCD collaborations have found stable NN- 
dibaryons and dibaryons containing s-quarks, but quark 
masses in their calculations are higher than the physi¬ 
cal values, see, e.g.^di. Chiral extrapolations of these 
results to the physical point gave, however, evidences 
against the existence of such dibaryons, see, e.g.?^. 
These calculations deal with ground states and say noth¬ 
ing about unstable states corresponding to a possibility 
of two-baryon fusion into 6-quark bag with a value of 
mass larger than a sum of masses of the initial baryons. 
Recent progress in excited baryon spectroscopy is de¬ 
picted i n^^d^ . Corresponding results based on nonphysi¬ 
cal quark masses too cover only one-baryon states so far 
and are in a poor agreement with experimental N and A 
excitation spectra. The first excited state in two-nucleon 
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system was found in lattice QCD ii>i^ but with a heavy 
quark mass corresponding to = 0.8 GeV. Therefore, 
predicting quasi-bound states of a multibaryon systems 
remains a difficult challenge in lattice QCD till now. 

Another important question: what is the reason that 
so few signs of dibaryons currently exist in spite of their 
search in the network of partial-wave analysis? The most 
likely answer, as we see it, is still an unsatisfactory preci¬ 
sion of PWA. Indeed, incorporating the additional data 
of WASA-at-COSY Collaboration into the SAID analysis 
produces a pole in support of the resonance hypothesis^^. 

A trivial generalization of the method of a cumulative 
particle is to select events with several, n > 1, secondary 
particles, not necessarily containing a cumulative one, 
which accompany a dibaryon production. Such a group 
of the additional particles, e.g., pions, may take away an 
excess of excitation energy, which put the main obstacle 
in the way of dibaryon recognition. In particular, Yu.A. 
Troyan reported the registration of some dibaryons us¬ 
ing just this methodii^— It should be, however, noted 
that most of experimental searches of dibaryons carried 
out in the past must be exposed to requalification. Let 
us consider, for example, a paper by B.M.Abramov et 
al3^, which is cited sometimes as a convincing argument 
against one of Yu.A. Troyan’s experiments. Even gross 
inspection of that paper reveals the following grave short¬ 
comings. Firstly, no methods of a background substruc¬ 
tion have been used. The solid line in the main figure 
of the paper— is only an optimal approximation of the 
experimental invariant mass spectrum containing, in the 
general case, a sum of background and dibaryon con¬ 
tributions. Secondly, number of events and precision of 
measurements do not allow to obtain a mass spectrum 
resolution nearly I MeV, which is necessary to verify 
confidently the Troyan results. Thirdly, the conditions 
of ” deep cooling”, which was ensured in Troyan’s exper¬ 
iments, has not been fulfilled in^S (see^ for details). 
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All dibaryons reported i n^^d^ were observed in in¬ 
elastic N-N interactions with additional secondary pi¬ 
ons. The extra pions take away an excess of excita¬ 
tion energy - a process which is a some kind of anneal¬ 
ing. This may reconcile two opposite requirements im¬ 
posed simultaneously on the system: it must be strongly 
compressed to form a compound state and it must be 
cold enough, since highly excited levels are usually short¬ 
living and elusive. Two additional pions in final state 
were in WASA-at-COSY and CELSIUS/WASA Collab¬ 
orations experiments^ d^d^ . Therefore, we may suggest 
with high reliability that synthesis of new multibaryons, 
and particularly dibaryons, should succeed an observa¬ 
tion made also for synthesis of new transuranium ele¬ 
ments - the system must be as much cold as possible to 
be observable readily. 
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For any vector n an identity (Sat • n)(SD • n) = 
^ ((S • n)^ — (Sjv • n)^ — (Sd ■ n)^) holds true. The term 
(Sjv • n)^ = 1/4 in the parentheses preserves S, the term 
(S-n)^ commutes with S^, but does not with S. This means 
that it conserves an absolute value of the total spin and 
breaks its direction. The term (So • n)^ does not maintain 
a direction of So and therefore a direction ofS = So + Sjv 
or even an absolute value of the total spin. 

This mathematical trick is described iir^ in section devoted 
to spin-orbit interaction. 

Hereafter we retain the overline as notation for averaging 
over Lorentz-invariant phase space. 

Using Cauchy-Bunyakovsky-Schwarz inequality, it may be 
proven that modulus of the interference terms defined by 
is in the general case greater than the true one. There¬ 
fore the role of interference is overestimated in our calcu¬ 
lations. Thus, we create an optimum for explanation of ex¬ 
perimental data by interference between different nucleon 
excitations, as far as the resonances have too large widths 
to explain cross-section by themselves. 



